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Abstract 

If the gravitino is sufficiently light and stable it will behave as an effective 
massless neutrino species at the time of nucleosynthesis. Depending on the 
temperature at which it decouples from the thermal bath in the early universe, 
the gravitino mass will be bounded by the primordial "^He abundance. Assum- 
ing a conservative estimate that the number of neutrino families, N^j < 3.6, 
superlight gravitinos with a mass 771.3/2 ~ lO^^eV are ruled out. This bound 
is weaker than previous estimates because the Goldstino annihilation cross 
section was overestimated. 
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I. INTRODUCTION 



If supersymmetry is to describe the real world then it must be broken at some scale A 
and give rise to supersymmetry breaking mass splittings not greater than m ~ 0{TeV). 
The signal for local supersymmetry breaking occurs when the superpartner of the graviton, 
the gravitino, receives a mass via the superHiggs mechanism. There is an upper bound of 
0{TeV) on the gravitino mass which comes from requiring that supersymmetry provides 
the solution to the naturalness problem. However, the gravitino mass is not necessarily 
constrained to be O(TeV) because the gravitino only interacts gravitationally and so is 
weakly coupled (compared to the gauge forces). This is not the case for the superpartners 
of the particles in the standard model which experience gauge forces. Their soft masses are 
constrained from experiment to be in the range (9(10^ — 10^)GeV. 

The relationship between the gravitino mass and the soft masses depends on the messen- 
ger sector. The messenger sector is responsible for communicating the spontaneous break- 
down of supersymmetry from a hidden sector, where the supersymmetry breaking dynamics 
occurs, to the observable sector. Normally in hidden sector scenarios of spontaneous super- 
symmetry breaking, the messenger sector is gravitational and gives rise to a gravitino mass 
m3/2 ~ m. A light gravitino is precluded because it would also mean light gauginos, sleptons 
and squarks. This typically assumes generic choices of the Kahler potential, Q{z,z*) and 
the gauge kinetic function fah{z) where z is the hidden sector chiral superfield. However, in 
no-scale supergravity models particular choices of Q{z,z*) and fah{.z) make it possible for 
the gravitino mass 7713/2 <^ m 

A hierarchy between 777.3/2 and fh is also possible if the messenger sector responsible for 
communicating supersymmetry breaking to the visible sector is not gravitational. Recently 
there has been a renewed interest in dynamical supersymmetry breaking models where 
gauge forces communicate the breakdown of supersymmetry to the squarks, sleptons and 
gauginos. In these models the soft masses {fh) are generated by radiative corrections and 
are proportional to the gauge couplings squared. The gauge sparticles will still receive soft 
masses from gravitational interactions, but these contributions are much smaller than the 
contributions from the gauge-mediated messenger sector. Since the gravitinos do not carry 
gauge quantum numbers they do not receive any mass from the gauge messenger sector 
but still receive mass from gravitational interactions. Consequently, in these models the 
gravitino is naturally light (-C (9(TeV)) and becomes the lightest supersymmetric particle 
(LSP). This can have interesting cosmological implications, such as contributing to the dark 
matter component of the universe |^. 

Given the possibility of non-minimal kinetic terms or a non-gravitational messenger 
sector one can ask how light can the gravitino mass be? A naive estimate of this mass 
comes from assuming that the supersymmetry breaking scale A ~ lO^GeV, which gives 
7773/2 ~ /Mpi ~ 10~^eV where Mpi is the Planck mass. If the gravitino were this light 
then it would essentially behave as a massless neutrino species during nucleosynthesis. This is 
because at energies E ^ 7713/2 the longitudinal component of the gravitino dominates during 
interactions, which is just a statement of the Equivalence Theorem [0. The longitudinal 
components of the gravitino are the helicity ±1/2 modes which come from absorbing the spin- 
1/2 Goldstone fermion (or Goldstino) during the spontaneous breakdown of supersymmetry. 
The Goldstino component of the gravitino couples with a strength proportional to K/7773/2, 
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where k = 1/Mpi. If the gravitino mass is hght this couphng can be much stronger than 
naively expected from the gravitational force. This means that as the temperature of the 
universe cools, the Goldstino will remain in thermal contact with the heat bath longer 
because of its enhanced coupling. Since the expansion rate of the universe depends on 
the number of degrees of freedom, the universe will expand faster with Goldstinos present 
at temperatures T ^ (9(100 MeV) and cause the neutrinos to decouple earlier. This will 
affect the production of neutrons (via the charged current weak interactions) by causing 
more neutrons to survive and ultimately increases the primordial ^He abundance. Thus the 
observed primordial ^He abundance sets a lower bound on the gravitino mass, which in turn 
is related to how late the Goldstino decouples from the thermal bath. 

A lower bound on the gravitino mass based on the nucleosynthesis argument was origi- 
nally discussed by Fayet ||^ , who considered the effect of the Goldstino together with a light 
photino on the effective number of massless neutrino species. The bound derived by Fayet 
(^3/2 ~ 10~^ eV) is not effective anymore because the photino is now believed to be of order 
the electroweak scale and the gravitino (if it is light enough) can only thermally interact 
with leptons and photons below the quark-hadron phase transition. More recently Moroi et 
al have also considered this bound and obtain an estimate m3/2 ~ 10~^eV. However the 
cross section used by them to recalculate the Fayet bound is only valid for T ^ m, In 
this work we will calculate the Goldstino decoupling temperature using the cross section for 
Goldstino annihilation in the limit m ^ T ^ ^3/2? where we assume that the supersym- 
metry breaking mechanism gives rise to a superlight gravitino. The gravitino mass bound 
will be shown to be much weaker than has previously been estimated. 

The outline of this paper is as follows. We begin in Sec. 2 with a discussion of the 
effective interaction Lagrangian for the Goldstino. This Lagrangian is used to calculate 
the Goldstino annihilation cross section into leptons and photons. In Sec. 3 the thermally 
averaged Goldstino annihilation rate is obtained. This will be used to calculate the Goldstino 
decoupling temperature in Sec. 4. A lower bound on the gravitino mass and the low energy 
supersymmetry breaking scale will be derived from the Goldstino contribution to the effective 
number of massless neutrino species. Final comments and the conclusion will be presented 
in Sec. 5. 



II. EFFECTIVE GOLDSTINO LAGRANGIAN 



The effective Lagrangian for the gravitino is obtained from the N=l supergravity theory 
0], irrespective of the mechanism employed in the messenger sector to communicate super- 
symmetry breaking to the observable sector. Denoting the chiral matter supermultiplets by 
(0, x) and the gauge vector supermultiplets by (A, A^) the supergravity Lagrangian for the 
interaction of the gravitino with the vector and chiral multiplets is given by 
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where a^'^ = f[7^,7'^], ^p denotes the gravitino and = iXi^Xi) is a Majorana spinor. 
The symmetric tensor h^^i, denotes the graviton where we have written the metric tensor 
as gpu{x) = rj^y + Kh^y{x). The hidden sector scalar components are defined as Sj = 
l/\/2{(f)j + 0*) and Pj = l/(-\/2z)(0j — (p*), and the constants Cj,dj satisfy 
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where nig is the gaugino mass, {fab) = {f)Sab and we have adopted the normahsing conven- 
tions used in Ref. [10]. In the Lagrangian (|l|) we have neglected all terms which are not 
relevant for the calculation of the gravitino annihilation cross section. 

The Lagrangian (|1|) may be used to calculate scattering processes involving all the helicity 
components of the massive gravitino. However, in the limit that the energy scale of the 
gravitino E ^ "^3/2? the longitudinal component of the gravitino will dominate and the 
gravitino effectively behaves as a spin-1/2 Goldstino 



' V3m3/2 



(3) 



where G denotes the Goldstino and the factor y| is a Clebsch- Gordon coefficient. This is 
just a statement of the Equivalence Theorem and is analogous to longitudinal W bosons in 
the standard electroweak model behaving as Nambu-Goldstone bosons in the high energy 
limit. If we substitute (H) into the supergravity Lagrangian (P and make use of the fact that 
for an on-shell gravitino 7'^\l/^ = then we obtain an effective Lagrangian for the interaction 
of the Goldstino with gauge and matter multiplets. 

The annihilation channels that will interest us are Goldstino scattering to leptons and 
photons. We will assume that all other particles have decoupled and are not important for the 
calculation of the Goldstino decoupling temperature. Consider first Goldstino annihilation 
into photons. There are three different contributions to this amplitude which are depicted in 
Fig. |l|. The first contribution is due to the t and w-channel exchange of the photino, where 
for simplicity we ignore neutralino mixing. In addition there are s-channel annihilation 
diagrams which result from graviton exchange and hidden sector scalar exchange. These 
diagrams are important for the cancellation of the leading order energy dependence in the 
Equivalence Theorem limit The Lagrangian relevant for these processes is given by 



^-/5(^MP^- + ^M<x^-P) _ p^^ pap 
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FIG. 1. Feynman diagrams for the scattering process GG — > 77. 
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Assuming that ^ ?Ti3/2 the Goldstino annihilation cross section into photons is calculated 
using the above Lagrangian (HI) to be 
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where x = m? / s. This agrees with the result quoted in Ref . [O . In the limit that ^/s ^ m. 



the cross section is in agreement with the Equivalence Theorem, namely 

(t{GG 77) 



1 4 4 
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where / generically denotes the photino "Yukawa" coupling. On the other hand in the 
intermediate limit ^ y/s the cross section (^) becomes 
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FIG. 2. Comparison of the various limits for the cross section a{GG 77) where only the 
quantity in square brackets (= ax) in Eq.^ is plotted. The solid curve represents the exact 
expression for ax, while the dot- dashed line depicts ax in the Equivalence Theorem limit ^ m) 
and the dashed line represents ax in the intermediate limit (rh ^ \^s). 



which has a different dependence on the energy scale y/s than in (H). This will be the limit 
that is needed to accurately calculate the Goldstino decoupling temperature. In Fig. ^ these 
various limits are compared with the exact result. It is clear from the figure that the cross 
section is a poor approximation in the limit m^y ^ ^/s. 

Notice also that the cross section (^) is proportional to E^. This is due to the non- 
renormalisibility of the supergravity Lagrangian and can lead to unitarity violation above a 
critical energy [|T^], 

-m3/2 



Ecr = V2887r— ^. (8) 



nrrig 



If the critical energy is taken to be the Planck scale, then it is difficult to reconcile a light 
gravitino mass with (H). However it was pointed out in Ref. [|12| that one can interpret 



the critical energy as the scale of new physics at which the gaugino mass, rrig becomes 
effective. This scale need not necessarily be anywhere near the Planck scale. For a gravitino 
mass m3/2 ~ (9(10~^eV) the critical energy turns out to be Ecr ~ 3TeV. This unitarity 



6 



FIG. 3. Feynman diagrams for the scattering process GG — > //. 

limit is safely above the energy scales that we will be concerned with, which are below the 
quark-hadron phase transition (~ (9(100 Me V)). 

The Goldstino annihilation amplitude into fermions is similarly obtained from the fol- 
lowing Lagrangian: 

C = t^^^id.Ghd'^h - hd.Gd^fl + fRd.Gd^a - d.Gfnd^fn) 
V o ms/2 

+ \^H>l^fd^Gl,rd,G + i-\f{Yd^ + l^d^)f - (^^/V + ^W)/] V 
b 777,3/2 4 

- - ^— {ri^^'r^P^ + r^^^r^P^" - r/^^r/^'') dM d^G h^p 

6 7773/2 

+ ^ ^ e"""'^" 75 [7<., dvG d^h^p 

b 7773/2 

- ^ E d^Irfd.P, -\^Y.d^ d,GYl5d^GdM (9) 
^ i ""^3/2 i 

where / denotes the fermion Dirac spinor. Note that due to gauge invariance there is no 
fermion coupling to the scalar Si. The relevant Feynman diagrams are depicted in Fig. |^. 
Besides the sfermion, graviton and hidden sector scalar exchange diagrams, there is also a 
four- Fermi interaction term due to the nonrenormalisability of the supergravity Lagrangian. 
Assuming that the fermions are relativistic the cross section is calculated to be 



a(GG ff) = ^m$s 

^ IOStT 771^/2 ^ 



1 + ^ — - i2x^ log 1 + 

I + X X 



(10) 
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FIG. 4. Comparison of the various limits for the cross section a{GG 
quantity in square brackets 



f f) where only the 
The solid curve represents the exact 



ax) in Eq.(lC) is plotted, 
expression for a^, while the dot- dashed line depicts ax in the Equivalence Theorem limit (a/s 3> m 
and the dashed line represents ax in the intermediate limit (rh ^ \^s). 



where x = m'^j/s and the contribution from the pseudoscalar P, diagram turns out to be 

zero. In the limit that ^/s nij one can easily check that (|1^) agrees with the Equivalence 
Theorem. If instead one considers the intermediate limit ^ i/i then the annihilation 
cross section (p!0| ) becomes 

aiGG ff) = — (11) 

where again the energy dependence is different than that in the high energy limit ^/s nij. 
The various limiting behaviours are shown in Fig. ^. One can clearly see in the figure that 
the Equivalence Theorem limit is not a good approximation when ^/s <^ m and infact the 
cross section is considerably smaller. 

To obtain the cross section for Goldstino decay into massless neutrinos one simply ne- 
glects the right handed components in the Lagrangian (^. In this case one finds 

a{GG ^ vlH) = \cr{GG ^ //). (12) 

The total Goldstino annihilation cross section is given by the sum of (^) and (0) and will 
be used to calculate the thermally averaged annihilation rate in the early universe. 
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III. THERMALLY AVERAGED ANNIHILATION RATE 



The annihilation cross sections obtained in the previous section were calculated at a 
temperature T = 0. In the early universe we need to average over the statistical distributions 
of the colliding particles in the thermal heat bath. The thermally averaged annihilation cross 
section times velocity for the scattering process 1 + 2 ^ F is given by |13| 



^^^^^^^^ = jdnVdnl^ ^^^^ 



where 



dnf = f{E,,t)g,^^ (14) 
and f{Ei, t) is the statistical distribution function and gi is the number of internal degrees of 



freedom for the particle species i. The factor vm01 = y {pi ■ V^)^ ~ '"^i^2/(-^i-^2) is known as 
the M0ller velocity (see Ref. [|1^) and a is the sum over all possible annihilation channels of 
particles 1 and 2. We have neglected the Pauli blocking factors 1 — / for final state particles 
and antiparticles in (p!3D, which typically amounts to a 10% correction in the determination 



of the decoupling temperature [14 



In the case of Goldstino annihilation we are assuming that E ^ ^3/2 so that the 
Goldstinos will be relativistic. Note that for Goldstinos which obey Fermi-Dirac statistics, 
f(E,t) = l/^e^f^ + 1) the equilibrium number density is 
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= ^C(3)r^ (15) 

where ({3) ^ 1.202. In the early universe we will be interested in the role of the Goldstino 
at temperatures T ~ 0(100 MeV). At these times <C fh and infact x = 10^ for fh ~ 
(9(100 GeV). If we parametrise the annihilation cross section as a = J2i (^iS^" where s is the 
Mandelstam variable and a, is a constant, then the thermally averaged annihilation cross 
section times velocity for the Goldstino in the limit T <^ m is given by 

(-«..) = 5^ (16) 

where /(n^) = /q°° dy ^5-^, ni are integers and the sum is over all annihilation channels. 

Using the expressions (^ and (0) for the Goldstino annihilation cross sections derived 
in the previous section we obtain 

= a{GG ^ 77) + E <GG ^ //) = a,,s^ + ^ a^js^ (17) 
/ / 

where we have summed over all possible fermion pairs in the final state. Assuming a photino 
mass ~ 0(100 GeV), the dominant part of the total cross section actually comes from 
the Goldstino annihilation into photons so that the thermally averaged cross section times 
velocity is approximately 
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{aAVM0i) ^ 1800^<T^^T^ (18) 



where C(5) — 1.037. The average Goldstino annihilation rate which is defined to be Ta = 
n^{aAVM0i) is then given by 

where a^^ is determined from Eq. (|^). This rate will be used to calculate the Goldstino 
decoupling temperature. 



IV. GOLDSTINO DECOUPLING TEMPERATURE 

During the radiation dominated era of the universe, the energy density which is domi- 
nated by relativistic particles, can be expressed in terms of the photon energy density p-yiT) 
as 

piT) = lg,iT)p,iT) = g,iT)^T\ (20) 

where the effective number of relativistic degrees of freedom of the bosons {B) and fermions 
(F) present is 

^.m^Es,{^)\l^,,{'f)\ (21) 



In Eq.(^ID gBigp) is the number of internal degrees of freedom for each boson (fermion) and 
Tb,f represent the possibility of the decoupled particles having a temperature which differs 
from the photon temperature, T. The Hubble expansion rate during this era will be 



The Goldstinos thermally decouple from the heat bath when their annihilation rate Ta ^ H. 
Using ([T9| ) and (p2D one finds that the Goldstino decoupling temperature is 
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This equation shows that as the gravitino mass becomes lighter, it causes the Goldstino 
to decouple later in the evolution of the universe. The later the Goldstino decouples, the 
greater the possibility of the Goldstino interfering with nucleosynthesis. Clearly there exists 
a lower bound on the gravitino mass for which the predictions of nucleosynthesis are not 
affected. 

During nucleosynthesis the energy density of new massless particles, i is equivalent to an 
effective number AiVj^ of additional doublet neutrinos: 
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(24) 



where = 8/7 for bosons, /i? = 1 for fermions and Qi is the number of internal degrees 



of freedom of the particle species i [1^. At neutrino decoupling the only known particles 
which can contribute to Qp are 7, e^, z/g, Ve-, ^n-, ^r, and v-r- This means that the effective 
number of degrees of freedom at neutrino decoupling are, using Eg. 

7 43 
^,(T,) = 2 + -(4 + 3x2) = -. (25) 

Assuming the conservative estimate that < 0.6 |16|, places a lower bound on the 

Goldstino decoupling temperature, T^. For example, if Td = T,^ — 0{MeV) then according 
to (p^), ANi, = 1 and the Goldstino would behave like an extra neutrino family. This would 
mean that Tjj > T^. A stronger bound can be obtained by supposing that the Goldstino 
decouples during the temperature range < < T^, where is the muon decoupling 
temperature. When the Goldstino decouples the effective number of degrees of freedom 
would be 

7 2^ 
^,(T,,) = 2 + -(4 + 3x2 + 2) = -. (26) 

Using Eq.(^) this will contribute an amount /S,Ny = 0.82, which means that the Goldstino 
decoupling temperature > ~ (9(100MeV). Imposing this condition on the expression 
(^) leads to the lower bound on the gravitino mass 

This bound is much weaker than that quoted by Moroi et al |^ because the Goldstino 
annihilation cross section is not as large as assumed by those authors. 

The mass bound (^) may be strengthened slightly by supposing that muons are also in 
thermal equilibrium when the Goldstino decouples. In this case one obtains gpiTr)) = 16 
and ANi, = 0.59. Assuming that the primordial '^He abundance rules this contribution out 
as well, causes the the lower bound on the Goldstino decoupling temperature to increase 



up to the pion mass T ~ m,r and the lower bound (^Tf) to increase slightly. Clearly, the 
higher the Goldstino decoupling temperature becomes the less it will contribute to N^, as 
many more particles and resonances contribute to gp{To)- In order to significantly increase 
the bound in future a more accurate estimate of AA^;^ would be needed. 

Of course it is likely that there are other beyond the standard model particles which could 
contribute to AN^. In this case it may be more difficult to accommodate the Goldstino, 



which would cause the bound (p7|) to increase further. This could happen for example if 



neutrinos have Dirac masses. It is known that neutrino Dirac masses can contribute a large 



amount to A A"^ . Accounting for the Goldstino and neutrino masses could place further 



constraints on the mass parameters. 
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V. CONCLUSION 



If a spontaneously broken N=l supergravity theory produces a sufficiently light gravitino 
then it will have consequences during the nucleosynthesis era of the early universe. This 
may be possible in no-scale supergravity theories or in theories of low energy dynamical 
supersymmetry breaking with a gauge-mediated messenger sector. The gravitino will obtain 
a mass via the superHiggs effect by absorbing the Goldstino. The enhanced coupling of the 
Goldstino causes it to interact more strongly with chiral and vector supermultiplets, which 
means that the Goldstino can decouple just prior to the nucleosynthesis era. 

The primordial ^He abundance critically depends on the number of massless neutrino 
families. If we require that the Goldstino not contribute significantly to the number of 
massless neutrino families, a lower bound on the gravitino mass can be obtained. Previous 
lower bounds have ranged from 10"^ — 10~^eV By calculating the Goldstino annihi- 

lation cross section into leptons and photons in the limit a/s <C m we were able to show 
that this bound is considerably weaker than previous estimates. If N^^ < 3.6 then typically 
^^3/2 ~ lO^^eV for ni:^ ~ (9(100 GeV). This bound complements previous gravitino mass 
bounds derived from collider experiments |T8| , p!9| and astrophysics . 

In addition a bound on the supersymmetry breaking scale, A can also be obtained for 
various scenarios of supersymmetry breaking dynamics. Assuming that m3/2 — A^/M^;, the 
bound on the gravitino mass implies that the scale of supersymmetry breaking, A ^ lOOGeV. 
In particular this would set a lower bound (100 GeV) on the scale of the supercolour sector 
in the recent gauge mediated models. These bounds are right at the forefront of existing 
collider energy scales. 
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